Magnetic Impurities in Mott-Hubbard Antiferromagnets 
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A formalism is developed to treat magnetic impurities in a Mott-Hubbard antiferromagnetic insu- 
lator within a representation involving multiple orbitals per site. Impurity scattering of magnons is 
found to be strong when the number of orbitals N' on impurity sites is different from the number 
N on host sites, leading to strong magnon damping and singular correction to low-energy magnon 
modes in two dimensions. The impurity-scattering-induced softening of magnon modes leads to 
enhancement in thermal excitation of magnons, and hence to a lowering of the Neel temperature 
in layered or three dimensional systems. Weak impurity scattering of magnons is obtained in the 
case N' = N, where the impurity is represented in terms of modified hopping strength, and a 
momentum-independent multiplicative renormalization of magnon energies is obtained, with the 
relative magnon damping decreasing to q 2 for long-wavelength modes. Split-off magnon modes 
are obtained when the impurity-host coupling is stronger, and implications are discussed for two- 
magnon Raman scattering. The mapping between antiferromagnets and superconductors is utilized 
to contrast formation of impurity-induced states. 
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I. INTRODUCTION 

While the problem of static impurities in antiferromag- 
netic insulators is more than twenty five years old, Q it 
has attracted renewed attention after the discovery of 
high-T c cuprate superconductors, || since their parent 
compounds are antiferromagnetic insulators. From the 
very early days of high-T c superconductivity a number 
of doping studies have been done with various static im- 
purities — both magnetic, || and nonmagnetic [||-|| - 
replacing copper from the Cu-0 planes as in La2Cu04. 
Susceptibility measurements in high-T c cuprates doped 
with magnetic impurities like Fe, Ni, Co give evidence of 
local-moment formation, Q which is expected to be in- 
trinsically associated with the magnetic impurities. This 
is unlike the case of nonmagnetic impurities such as Zn, 
Al, Ga etc. which, despite being intrinsically nonmag- 
netic, give rise to local moments in the copper-oxide 
planes when doped in cuprate antiferromagnets. This 
was inferred earlier from the Curie- Weiss behavior of the 
magnetic susceptibility, jJJ^] and has been recently con- 
firmed in the Y-NMR studies of doped 1-2-3 systems as 
seen in the progressively increasing lincwidth of the Y- 
NMR signal with decreasing temperature. J|,|| Xiao et 
al. have also ascertained the spin states of different mag- 
netic dopants from the observed local moments, || and 
find, for example, that Fe is in a spin-| state, whereas Ni 
is in a spin-1 state. They also find a correlation between 
T c reduction and size of the local moment, consistent 
with the magnetic pair breaking mechanism. 

Although theoretically the problem of magnetic im- 
purities in an antiferromagnet has been studied recently 
within the Heisenberg representation of localized spins, 
Q no such comprehensive study exists within the Mott- 
Hubbard model, which provides a good description of the 



3d holes in the Cu-0 planes of high-T c superconductors. 
Recently the problem of nonmagnetic impurities in the 
Mott-Hubbard antiferromagnet was addressed and de- 
fect states, local-moment formation, impurity-scattering 
of magnons, and finite-temperature magnetic dynamics 
in layered systems were studied. Other recent 

works on static vacancies in antiferromagnets include ex- 
act diagonalization studies with Heisenberg model, 
linear spin wave theory, Jl3| and exact diagonalization 
of strongly correlated small clusters. 0] While non- 
magnetic impurities can be simply represented by spin- 
independent impurity potential, the situation is more 
complex for magnetic impurities. In this paper we in- 
troduce several representations to treat magnetic impu- 
rities in different situations. A simple extension to spin- 
dependent impurity potential is followed by a more so- 
phisticated approach involving a generalized A/"-orbital 
Hubbard model with multiple orbitals per site. Broadly 
there are two distinct classes depending on whether the 
number of orbitals N at the impurity site is the same 
as or different from the number of orbitals N at the host 
sites. In the case N = N the magnetic impurity is rep- 
resented through a modified hopping strength if between 
the impurity orbitals and the neighboring host orbitals. 
In the strong-correlation limit (U 3> t) wherein the Mott- 
Hubbard AF with N orbitals per site maps to the spin 
S = N/2 quantum Heisenberg AF, the modified hop- 
ping strength translates into modified exchange coupling 
J' = 4t' 2 /U between the impurity spin and the neigh- 
boring host spins. This describes the situation where, in 
spin language, the impurity spin S' is equal to the host 
spin S, but is coupled to its neighbors with a different ex- 
change interaction J'. Similarly the case N ^ N with no 
modification in hopping strength or Hubbard interaction 
energy corresponds to the situation where the impurity 
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spin is different from the host spins (S' ^ S). 



II. SINGLE-ORBITAL MAGNETIC IMPURITY 

In this section we consider a single-orbital magnetic im- 
purity embedded in an AF host which is described by the 
Hubbard model with one orbital per site with exactly half 
filling. For concreteness we consider the square lattice, 
generalization to other bipartite lattices being straight- 
forward. The host Hamiltonian is 



(1) 



where t is the nearest-neighbor (NN) hopping strength 
and U the on-site Coulomb repulsion. The AF state and 
its associated features such as sublattice magnetization, 
magnon energies, quantum corrections etc. have been 
studied earlier in detail. 15 We model the single-orbital 
impurity in terms of locally modified hopping term t' 
between the impurity orbital and its NN host orbitals. 
The Hamiltonian with such an impurity on site / can 
be written as below, where the sum is over all nearest 
neighbors J of the impurity sites /, and St = t' — t is the 
hopping perturbation around the impurity site, 



H = Hn 
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We start with the perturbative method where the 



impurity-induced perturbation [6x°] = [x°] — [Xh 



host- 



to 



the zeroth-order, antiparallel-spin, particle-hole propa- 
gator is obtained in powers of St/t, and resulting cor- 
rections to its eigenvalues then yield the renormaliza- 
tion in magnon energies. |16| Diagrammatic contribu- 
tions to [Sx°] to first order in St, and their evaluation in 
the strong-correlation limit have been discussed earlier in 
context of the hopping disorder problem. |l(J We obtain 
for the diagonal, off-diagonal, and nearest-neighbor diag- 
onal terms, expressed in units of — t 2 / A 3 for convenience, 

[Sx°]n = f j; [Sx°}ij = [Sx°hi = [Sx°hj = \j 

(3) 

where z = 4 is the coordination number for the square 
lattice, 2A » U is the Hubbard gap, and only terms 
upto order (t 2 /A 3 ) have been retained, appropriate to 
the strong-correlation limit. We notice that the sum of 
all matrix elements diagonal in sublattice basis, [#x ]jj + 
[^X°] J J is precisely equal to the sum of off-diagonal ma- 
trix elements [<5x°W + [^X°]j/- An immediate conse- 
quence of this correlation is that the Goldstone mode is 
preserved and that generally the effective scattering of 
low-energy, long-wavelength magnon modes is weak. 

If the impurity is on an A-sublattice site, then for the 
first-order correction we obtain after summing over near- 
est neighbor terms, 



S\W = {q\[S X }\q) = 

a 2 [S X °]n + aftz lq [S X %j + ftaz lq [S X \ u + ft 2 z[S X °]jj (4) 

where a and ft are the magnon amplitudes on A and 
B sublattices respectively, and j q = (cosq x + cosq y )/2. 
An identical result is obtained when the impurity is 
on a B-sublattice site, because in this case a and ft 
are simply exchanged in the above equation, and since 
[^X°]i7 = z [^X°}jJj this expression is symmetric under 



where u>„ 



exchange of a and ft. Using a = ^Jjj(l— cj°) and 



ft 



1 — 7? is the host 



magnon energy in units of 2 J for the momentum-g mode, 
we obtain after summing over contributions from all im- 
purities 



(5) 



where x is the total impurity concentration, and impuri- 
ties are assumed to be evenly distributed between the two 
sublattices. The renormalized magnon energy, given by 
the pole in the magnon propagator, is now obtained from 

the solution of the equation 1 — ^Jlu 2 + ^ 2 + 5X^ = 0, 

and upto first order in the effective impurity strength 
xSt/t we obtain 



St 

u^yl + xz— 



,o 



(6) 



This result agrees exactly with the calculations |9| on 
the Heisenberg model in that there are no singular cor- 
rections to the magnon energy in the case 5" = 5 1 , and 
the correction is proportional to xSt/t = (l/2)xSJ/J. 

Turning now to the magnon-energy renormalization of 
the localized, high-energy modes with energy near 2J, 
which correspond to local spin deviation, we have a = 0, 
ft = 1, so that <5A (1) = \z{5t/t). This implies that the 
magnon energy gets shifted from 2 J to 
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In this case the impurity concentration does not enter the 
magnon-energy renormalization, rather it has a bearing 
on the spectral weight of these high-energy modes. Thus 
for St positive, the magnon spectrum goes up by energy 
z J St/t. This increase is expected from the simple picture 
of these high-energy modes corresponding to local spin 
deviations. The energy cost of making a spin deviation 
on the impurity site is zJ'/2, where J'/2 is the bond 
strength. With t' = t(l + St/t), to first order in St/t we 
have Ae = z(J' - J)/2 = z J St/t. 

The exact-eigenstates analysis also shows that pre- 
cisely one magnon state at the upper end of the spec- 
trum is split off from the magnon energy band. These 
split-off modes are strongly localized around the impu- 
rity sites, and hence correspond to local spin deviations. 
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Furthermore, for different values of the impurity hop- 
ping t' jt it is seen from the magnon spectrum that the 
energy separation of the split-off state from the upper 
end of the spectrum increases roughly in proportion to 
St, as obtained in the perturbative analysis. This exact- 
eigenstates approach for obtaining magnon energies and 
wavefunctions from the fermionic eigensolutions in the 
self-consistent AF state has been described earlier. jlTj 

III. SPIN-DEPENDENT IMPURITY POTENTIAL 

Nonmagnetic impurities in the Mott-Hubbard AF were 
modelled earlier via a spin-independent impurity po- 
tential term, and as a natural extension we therefore 
consider the following spin-dependent impurity term for 
magnetic impurities, 

fl££ = £*/[-"3V]*/, (8) 

where tyj = (d/f djj.)- A spin-independent impurity 
potential eo can be included for generality, however, we 
shall consider the limit V >> eo, so that the potential for 
spin a fermion is V a ~ —crV. We choose V to be positive 
for impurities on the A-sublattice sites, so that V\ is very 
low and V± is very high. The sign of V is reversed for im- 
purities on B-sublattice sites. This choice of potential en- 
sures that the magnetization on the impurity sites follows 
the host AF ordering. Such a spin-dependent impurity 
potential can arise from a coupling — cr.Si mp between the 
itinerant fermion spin a and the static magnetic impurity 
spin /f>i mp , resulting from a strong Hubbard interaction. 
Since experiments on high-T c cuprates show the impurity 
spin to be antiferromagnetically coupled with the host 
spins, ^ we take the local field direction to be along the 
local magnetization direction (£). The low potential (for 
spin-up) is justified in view of the fact that the ioniza- 
tion energy for both Fe +3 and Ni +2 i.e., the fourth and 
the third ionization energies respectively for Fe and Ni 
are much higher than the third ionization energy for Cu. 
Whereas the ionization energies for Fe +3 and Ni +2 are 
54.8 eV and 35.17 eV respectively, the ionization energy 
for Cu+ 2 is 20.2 eV. 

We now examine formation of impurity-induced states 
due to this spin-dependent impurity potential. Within 
the T-matrix analysis, used earlier for nonmagnetic im- 
purities, []lof energies of impurity-induced states are ob- 
tained from solutions of gfj(oj) = 1/V a . For large \V\/U 
these impurity states are formed at energies ~ —aV for 
the two spins, and are essentially site localized and there- 
fore decoupled from the system. Thus, for the magnetic- 
impurity case when the impurity spin is antiferromag- 
netically coupled to the neighboring host spins, a sig- 
nificant difference from the nonmagnetic-impurity case 
is that there are no defect states formed in the Hubbard 
gap. Rather only impurity states are formed, far removed 
in energy from the Hubbard bands. 



Within the above representation of magnetic impuri- 
ties in terms of spin-dependent impurity potential, the 
fermion number is unchanged, unlike the case of nonmag- 
netic impurities where one fermion is removed for every 
added impurity. Hence the impurity sites do not quite act 
as spin vacancies. Nonetheless, the presence of a impu- 
rity potential term which breaks time-reversal symmetry 
leads to a partial decoupling of the impurity site from the 
host. This is most easily seen in the limit V — > oo where 
the local antiparallel-spin, particle-hole excitations are 
suppressed by the large energy difference 2V, leading to 
an absence of the u) term, and therefore to strong magnon 
scattering. Quite generally, the particle-hole energy dif- 
ference for antiparallel spins is modified by the spin- 
dependent impurity potential from 2 A to 2A + 2V, lead- 
ing to a modification in the w term. For spin-independent 
impurity potential the particle-hole energies are shifted 
equally, and hence it is the removal of a fermion from the 
impurity site that is crucial. As a result of this decou- 
pling of magnetic impurity sites, a qualitatively identi- 
cal impurity-induced perturbation is obtained, 
leading to similar results for magnon renormalization as 
for the nonmagnetic-impurity case, where singular cor- 
rections were obtained for low-energy magnon modes in 
two dimensions, jnj The strong impurity-scattering of 
magnons also introduces significant damping, with the 
ratio of magnon damping term to its energy being sim- 
ply proportional to the impurity concentration x for long- 
wavelength modes. 

IV. GENERALIZED HUBBARD-MODEL 
REPRESENTATION 

In order to represent higher-spin magnetic impurities, 
we now generalize to the situation with J\f and A/ 7 or- 
bitals on host and impurity sites respectively. An ap- 
propriate model for this case is the generalized Hubbard 
model with multiple orbitals per site. This model has 
been used earlier to study quantum corrections in the 
antiferromagnetic state in a spin-rotationally-symmetric 
formalism, where a systematic per turbative expansion in 
powers of 1/Af was developed. []15j We introduce a slight 
extension here in this model which makes it equivalent, in 
the strong correlation limit, to the spin-S* QHAF, where 
S = AF/2, The modification is to allow the NN hopping 
term to operate between all orbitals, whereas the hopping 
term considered earlier was diagonal in the orbital index. 
|l5| We therefore consider the following Hamiltonian for 
the AF host, 

H = ~t ( a La«ja/3+h.C.) 

<ij>aa@ 

+ Jf £(4Ta a *ta a i|/3 a il/3 + a h a a *w4lf3 a *lc*) (9) 

where a and (5 are the orbital indices which run from 
1 to A/", and the two Hubbard interaction terms are re- 
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spectively direct and exchange type interactions with re- 
spect to orbital indices. In the symmetric case when 
the two interaction strengths are identical, as consid- 
ered here, the system possesses spin-rotational symme- 
try. It has been shown earlier that in the symmetric 
case the two interaction terms can together be written 
as Hi nt = -{U/Af)J2i(Si-Si + n i)> where Si and 
are the total spin and charge density operators, respec- 
tively. Spin-rotational-symmetry is therefore inherent in 
this impurity representation as well. Furthermore, in the 
strong correlation limit, a strong Hund's coupling exists 
which energetically favors the maximum multiplicity case 
(S = Af/2) for the total spin operator Si. 

Magnetic impurities are represented by introducing 
TV' 7^ Af orbitals at the impurity sites. We first examine 
the transverse spin fluctuation propagator in the host 
AF state, x~ + (rt,r't') = ($\S-(rt)S + (r't')\$), where 

S( r t) — E Q Vd( r 0§ V ; a( r is the total spin operator. 
Again, at the RPA level the magnon propagator is given 
by — (C//A/")x°(w)], where [x°( w )] now involves 

orbital summations, with matrix elements given by, 



\x»]i 



Q/3 



(wOG^ia^'-w)- (10) 



Since each orbital is now connected via hopping to 
Af orbitals on the NN sites, the electronic spectral 
weights are correspondingly modified. For example, in 
the strong-correlation limit, the on-site majority and mi- 
nority spin densities in each orbital are now 1 — Aft 2 / A 2 
and Aft 2 / A 2 respectively. A straightforward extension 
of the earlier analysis in the strong-correlation limit |l5| ] 
leads to: 



1 



D t 2 



x »{ q ,u J )=N-l-N 2 -^ 



U~ 



1 



DJN 



lq 

'DJN' 



(ii) 



where J = 4t 2 /U as usual, and D is dimensionality 
of the hypercubic lattice. Since different orbitals on 
the same site are not directly coupled, the intrasite 
propagator is diagonal in orbital index, and therefore 
the leading order diagonal terms (the 1/U and the w 
term) are proportional to Af. However, the NN hop- 
ping operates between all orbitals, and therefore the off- 
diagonal term and the next-to-leading order piece (aris- 
ing from hopping) in the diagonal term are both propor- 
tional to Af 2 . The magnon energies are now given by: 

u q = DJN J I -j 2 = 2DJSjl-rf in terms of the 



spin S = Af/2. 

We now introduce a magnetic impurity in the system 
with Af' ^ Af orbitals at the impurity site /. The result- 
ing modification in the electronic spectral weights leads 
to the following changes in the [x°( w )]y - matrix elements 
for i,j in the vicinity of the impurity site I : 



DJAfJ 



u 



-AfAf' 



Df_l 
"2A 3 ! 



i 1 DJAf) 



■ Af(Af' - AT) 



D t 2 1 



(12) 



2 A 3 z 

Since now the local Hubbard interaction strength itself 
is not uniform but depends on the number of site orbitals, 
we have to multiply the [% ] matrix with the diagonal in- 
teraction matrix [U] containing elements U /Af for host 
sites and U /Af' for the impurity site. We therefore ex- 
amine the local matrix elements of the matrix product 
\Ux ]ij = MaXij f° r hJ m the vicinity of the impurity 
site. The impurity-induced perturbation in the matrix 
elements of the product [Ux°] are obtained as below: 

5[U X °{u)]u =0 
8[Ux q {oj)]ij =0 

Spx^U = -U(AT'-Af) D '' ' 



2 A 3 2 

S[Ux»]jj =-U(Af'-Af)^- z . 



(13) 



We now obtain the magnon-energy renormalization 
by perturbatively obtaining the impurity-induced cor- 
rection to the eigenvalues of the [Wx°(o;)] matrix. As 
discussed earlier, [[n] we treat 5[Ux°(u)] as the pertur- 
bation matrix, and determine corrections to eigenval- 
ues of [Xhost( u; )]- Evaluating the first-order correction 
((7|(5[Wx°(t<-')] \q) from the magnon eigenvector |g), and re- 
taining terms to first order only, we obtain: 



SAW 



U(Af' - Af) 



D t 2 



(14) 



2 A 3 DJ' 

As for the nonmagnetic impurity case, we obtain here 
a correction which is linear in energy, and this sig- 
nifies strong impurity scattering of magnons for long- 
wavelength, low-energy modes, leading to singular cor- 
rections in two dimensions and strong magnon damping 
from second-order scattering processes. |nj The scat- 
tering term is explicitly proportional to the difference 
(Af' — Af) between the number of orbitals on the impu- 
rity site and the host sites, which arises from the different 
dynamics of the impurity spin and the host spins. This 
generally implies that impurity scattering of magnons is 
strong when the impurity spin S' = Af' /2 is different from 
the host spin S — Af/2, in agreement with earlier studies 
within the Heisenberg model, || and the one-band model 
with nonmagnetic impurities where Af = 1 and Af' = 0. 

m 



V. CONCLUSIONS 

In conclusion, we have developed a formalism to treat 
magnetic impurities in a Mott-Hubbard antiferromag- 
netic insulator within a representation involving multi- 
ple orbitals per site. For the case Af' = Af, when the 
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impurity spin is identical to the host spin, the mag- 
netic impurity is represented by locally modified hop- 
ping strength, and we find that the effective scatter- 
ing of long-wavelength magnon modes is weak, lead- 
ing to momentum-independent multiplicative renormal- 
ization of magnon energies. For positive hopping per- 
turbation 5t we find localized, split-off magnon modes 
corresponding to local spin deviations at impurity sites. 
These split-off modes will be relevant in two-magnon Ra- 
man scattering which probes high energy magnetic exci- 
tations. In the other case A/"' ^ TV, when the impurity 
spin is different from the host spin, we obtain strong 
impurity scattering of magnon modes proportional to 
the difference (A/ 7 — TV) , leading to singular corrections 
in two dimensions and strong magnon damping. The 
impurity-scattering-induced softening of magnon modes 
implies enhancement in thermal excitation of magnons, 
and hence to a lowering of the Neel temperature in lay- 
ered or three dimensional systems. We also find that 
the process of putting additional impurity orbitals leads 
to enhanced impurity magnetization and localization of 
electronic states at the impurity, indicating partial de- 
coupling of the impurity site from the host. A unique 
feature of having multiple impurity orbitals is the pres- 
ence of exactly site-localized eigenstates in the electron 
spectrum which are completely antisymmetric between 
impurity orbitals. 

When the magnetic impurity is represented in terms 
of a spin-dependent impurity potential, we find that the 
breaking of time-reversal symmetry leads to a decoupling 
of the impurity site from the host, and strong magnon 
scattering similar to the case of spin vacancies is ob- 
tained. We also find that when the magnetic impurity 
spin is antiferromagnetically coupled to the neighboring 
host spins, only impurity states are formed, and there are 
no defect states formed within the Hubbard gap. The lo- 
cal moment associated with the magnetic impurity there- 
fore intrinsically arises from the spin-density difference 
at the impurity site. Using the well-known particle-hole 
transformation, The problem of magnetic impurities in 
an AF can be mapped to that of nonmagnetic impurities 
in a superconductor, which is characterized by absence 
of defect states within the superconducting gap and ro- 
bustness of superconducting gap. [|l9j Conversely, a non- 
magnetic impurity in a positive-t/ Hubard AF maps onto 
a magnetic impurity in a negative-C/ Hubbard supercon- 
ductor, and here defect states are formed within the gap 
in both cases. [fo]|2l) 
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